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Introduction. 

§ 1. In a former paper # I have considered the effect of gravitation as a factor tending 
towards instability, in the case of a spherical nebula of gas. The object of the present 
paper is to investigate the analogous problem in the case of a spherical planet, the 
planet being supposed composed of solid or fluid matter. The main question at issue 
is the following. 

§ 2. So long as gravitation is neglected there can be no doubt as to the stability of 
an elastic solid; any displacement increases the potential energy, and an unstressed 
configuration of equilibrium is therefore necessarily stable. But when gravitation is 
taken into account, the gravitational energy may be either increased or decreased by 
a displacement from equilibrium, and if a displacement can be found which effects a 
decrease in the gravitational potential energy of amount sufficient to outweigh the 
increase in the potential of the elastic forces, then the equilibrium configuration will 
be unstable. 

Now, in § 2 of the previous paper already referred to, it was shown that for any 
spherical body displacements can be found such that there is a decrease in the 
gravitational potential. This is sufficient to prove that a spherical configuration of 
equilibrium may be unstable. 

In the terminology of PoiNCAREf it appears that on any " linear series" of 
spherical configurations there may be " points of bifurcation." 

We must, therefore, attempt to settle the position of these points of bifurcation. 

In particular, it will be of interest to examine whether a sphere of the size and 
material of the earth may be regarded as being anywhere in the neighbourhood of a 
point of bifurcation. 

* " The Stability of a Spherical Nebula," < Phil. Trans.,' A, vol. 199, p. 1. 
t ' Acta. Math./ vol. 7, p. 259. 
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Prelimin ary Approximation. 

§ 3. A rough and very simple calculation will give an approximate answer to this 
latter question. 

Let a be the radius of a sphere, which will ultimately be taken to be the earth, 
M its mass, and p the mean density given by M = %irp Q a s . 

Let us use the elastic constants X, /a* and let X be the mean value of X. Since 
the sphere is supposed to be spherically symmetrical, X, ja, and p will be functions of 
the single co-ordinate r, the distance from the centre. Imagine X/X , /x/X , and p/p 
each expressed as functions of r/a, and let c l9 c 2i . . . be the coefficients which occur 
in these functions, these coefficients being mere numbers and independent of the 
system of units in which X, p, and a are measured. 

Imagine a linear series of equilibrium configurations obtained by varying any one 
of the quantities X , p , or a, while keeping the remaining two quantities and the 
coefficients c 3 , c 2 , . . . constants. The points of bifurcation on this series will occur 
when the varying parameter becomes equal to some definite function of the remaining 
quantities and of y, the gravitational constant. 

Hence, however the linear series are arrived at, the points of bifurcation will be 
iven by an equation of the form 



g 



J \Ys ^-o? Po? G&j £i> Cg, . . .j — u (1). 

Now the coefficients c 1? c 3 , . . . are mere numbers, and the only way in which y, 
X , p Q , and a can be combined so as to give a mere number is through the term 
yp 3 a 2 /X . Hence equation (1) can be expressed in the form 

f (Wo*** „ \ _ a (9\ 

/I , t^lj ^Oj ... I — — Vj ........ \£J l# 

We have seen that the spherical configuration must be unstable for some values of 
y, p , a, and X (e.g., it is always unstable for yp 2 a 2 /X = oo ), hence equation (2) must 

have at least one real root between yp 2 a 2 /X — ® anc ^ TPg^Ao ' =z °° • -^ e ^ the l° wes t 

root be 

y/) s a 2 /X = <j) . . . . . . . . . . (3), 

where <£ is a function of c 1? c 2 , . . . ; then a spherical configuration is stable so long as 
yp 3 a 2 /X < cj>, and becomes unstable as soon as yp 2 a*/\ > <£. 

The coefficients c l9 c 2 , . . . will, on the average, be comparable with unity, because 
X, p are referred to their mean values ; they are as likely (speaking somewhat loosely) 
to be above as to be below unity. Hence cf> itself will be comparable with unity, and 

* The notation is that of Love's ' Theory of Elasticity.' The m, n of Thomson and Tait are given 

A + fi = m, n = n. 
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it is not at present possible to say whether it is more likely to be greater or less 
than unity. 

§ 4. Now, in the case of the earth (Thomson and Tait, § 838), we have 

a = 640 X 10 6 centims., p = 5*5, 

and the value of y in C.G.S. units is known to be 

y = 648 X 10~ 10 . 

This gives for yp^a 2 the value 

y/) V=8.Xl0 u , 

whence it appears that for a sphere of the size and mass of the earth the spherical 
configuration will be unstable unless X has a value comparable with 8 X 10 11 . 

Now for steel (cf. Thomson and Tait, p. 435) the values of the elastic constants in 
absolute units are n = fx = 7 '7 X 10 11 , m = X + p = 16*0 X 10 11 , whence 
X = -8*3 X 10 11 . We therefore see that the critical values of the elastic constants in 
the case of the earth are comparable with those of steel. 

The foregoing calculation is, of course, very rough, but it shows that the critical 
values for the earth are at least in the neighbourhood of what must be supposed to be 
the actual values, so that we are driven to attempting a more accurate determination 
of these values. If the view of the present paper is sound, this approximate equality 
is more than a mere coincidence ; we shall see that it could have been predicted from 
our hypotheses of planetary evolution. 

We now attempt a rigorous mathematical investigation of certain problems which 
have a bearing upon the astronomical questions in hand. Those readers whose 
interest lies in the application of the results rather than in the processes by which 
they are obtained may be recommended to turn at once to § 22. 



The Stability of a Gravitating Elastic Solid. 
The Equations of Small Vibrations. 

§ 5. We shall begin by discussing the principal vibrations and the frequency 
equation of a spherically symmetrical solid. The case of a non-gravitating sphere 
has been fully discussed by Professor Lamb, # but the inclusion of the gravitational 
terms, as will be seen later, brings about a considerable complication in the analysis. 
The case of a gravitating but incompressible sphere has been considered by 
BROMWiCH,f but this has no bearing on the present problem, in which the whole 

* "On the Vibrations of an Elastic Sphere," 'Proc. Lond. Math. Soc.,' vol. 13, p. 189. 

t "On the Influence of Gravity on Elastic Waves, &c," £ Proc. Lond. Math. Soc./vol. 30, p. 98. 
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interest turns upon the compressibility. The solution which follows is, in its main 
points, very similar to that of Professor Lamb, so that I have not thought it 
necessary to give the steps of the argument in great detail. 

From the symmetry of the solid it follows that the elastic constants X, /x, and the 
density />, will be functions of the single co-ordinate r> the distance from the centre. 
Taking the centre as origin, we shall use rectangular co-ordinates, x, y, z, and shall 
suppose the solid to execute a small vibration, such that the displacement of the 
element initially at x, y> z has components, £ rj, £. The component of displacement 
along the radius will be denoted by u and the cubical dilatation by A, so that 

u = - (& + yy + £s), A = ~~-+ t+7° 

r x J ax ay dz 

§ 6. After displacement the density at x, y, z is 



or, since p is a function of r only, 



dp 

p _ Ap _ tt _ 



Hence the gravitational potential at x, y, z is changed by displacement from V 
into V — E, where E is the potential of the following distribution of matter :— 

(i.) A volume distribution of density 

. dp 

(ii.) A surface distribution of which the surface density is 

this being taken over every surface at which the density changes abruptly, the change 
being from p to p i in crossing the surface in the direction of r increasing. In 
particular this will occur at the outer surface of the solid, the value of p i in this case 
being zero. # 

§ 7. The potential at x, y, z after displacement being V — E, that at x + £, y + y), 

z + £ will be 

3V 3V 3V t 3 3 Y 

?&? ^ ^3? '•• . 

* In the investigations on gravitating spheres given in Thomson and Tait's < Natural Philosophy/ the 
course of procedure is tantamount to neglecting the volume distribution (4), and regarding E as the 
potential of a surface distribution (5) alone. For this reason the result obtained differs from that of the 
present paper. 
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Hence the force at x + L V + V> z + £ m the direction of x increasing, found by 
differentiating the foregoing expression with respect to £, is, neglecting squares of the 
displacements, 



3V t c 3 V 3-V 3"V 8E 

a, : + ^ a^' + v a"r 0// + ^ a* : & " a*" ••■•••• ( 6 )- 

Let us suppose that, in addition to its own gravitation, the sphere is acted upon 
by an external field of force of potential V , and let us, in the usual notation, denote 
the six stresses by P, Q, II, S, T, U. Then the equations of motion of the element 
at x + f ? y + y], % + £ i* 1 the displaced configuration are three of the form 

3? f — 3r j_ 3IT i 3T i / 3W _l £ 3 ' W j, 3 " W a. / S3W _ 3E \ 

p ¥ " a.; + a7 + & + p V"a7 + * a^ + v a* d y + * a7a* & / • * ( 7 '> 

- . aw 3 2 w 

in which W = V + V n , and all the terms such as -~ — , -tttt j • • • are evaluated at 
x, ?/, z, but p, P, Q . . . are calculated in the displaced configuration at x + f, y + ?), 

§ 8. Now r the only case in which we have any accurate knowledge as to the values 
of P, Q, R, S, T, U is when the whole strain is small, i.e., when W is small. In the 
case of the earth, V is not, in this sense, small. # The only way in which we can 
proceed with any certainty is, therefore, by taking V = — V, or W = 0. That is 
to say, we must artificially annul gravitation in the equilibrium configuration, so that 
this equilibrium configuration may be completely unstressed, and each element of 
matter be in its normal state. In this case it seems justifiable to suppose both the 
density and rigidity to be constant throughout the sphere, and, indeed, it is only 
with the help of this simplification that the equations become at all manageable. 

The vibrations of this system will be of two kinds. First there are "spherical" 
vibrations in which the displacement is purely radial at every point, so that the solid 
remains spherically symmetrical after displacement, and, secondly, there is the larger 
class of vibrations in which the displacement is not of this simple type, so that the 
displaced configuration is not one of spherical symmetry. 

We hope, by discussing the vibrations of this system, to obtain some insight into 
the corresponding vibrations of a natural non-homogeneous solid, say the earth. Now 
it is extremely doubtful whether the spherical vibrations of our artificial system have 
much in common with those of the natural system, but it will be seen later that this 
is of no importance. We shall not be in any way concerned with these vibrations. 
What we shall require is a knowledge of the unsymmetrical vibrations, and this, it is 
hoped, can be obtained with fair accuracy from a consideration of the corresponding 
vibrations in the artificial case. There must be some uncertainty even in the case of 
unsymmetrical vibrations, and, unfortunately, this seems to be inevitable ; our 



* Love, ' Elasticity/ I., p. 220, 
VOL. CCL— A. Y 
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artificial case appears to be the only case in which the equations can be solved by 
ordinary analysis. 

We now replace P, Q, R, S, T, U by their ordinarily assumed values, and equation 
(7), putting W = 0, takes the form 

P^ ' = ( X + /*)&T + ^ v2 f ~P^' * •••••• ( 8 )> 

and there are two similar equations for 77, £. 

The Principal Vibrations and Frequency Equations. 

§ 9. Differentiate these three equations of motion with respect to x, y, z and add ; 
then 

p ^ = (X + 2/x) V 2 A - pV 3 E (9). 

Now, from the definition of E, we have, in the case in which p is constant, 

V 2 E = ~~ 4.<rrpA . . . (10), 

and hence equation (9) becomes 

•72A 

p -— = (X + 2^) V 2 A + 4tt/) 3 A ....... (11). 

If we suppose A proportional to cos pt, this equation assumes the form 
(V s -f- ^) A = 0, where 

2 _ p(^ 3 + 4tT/)) / x 

fV — — — „ e «, 9 . . * . r V / 



There is, therefore, a particular solution of (11) of the form 

A = r~ i J n + i(Kr)8 /l (0, <f>) cos pt ....... (13), 

where S» (0, </>) is a surface harmonic of order n 9 and the general solution found by 
summation of solutions of this type is 

A = S2r"*J„ + i (kt) S„ (0, </>) (A cos jp* + B sin pt) . . . . (U), 

where the summation extends over all possible harmonics, and over all values of k. 

It will appear later that each term in this solution can be made to satisfy the 
boundary conditions, and, therefore, that each term represents a normal vibration. 

The vibrations may, therefore, be classified into vibrations of order 0, 1, 2, &c, the 
order being that of the harmonic which occurs in the expression for A. The vibrations 
of order n = are the spherical vibrations already referred to. 
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We shall assume this provisionally, in order to avoid the continual repetition of 
double summation, and now proceed to evaluate f, ??, £ and to form the boundary 
equations for the simple vibration given by equation (13). 

§ 10. From equation (8) it appears that the displacement ^ is given by 

p 3 pf + /xV 2 £ = - (X + /x) ^ + p -^ (15). 

The solution is 

where <£ is any solution of 

P*p<f> + /^V 2 (j& = — (A. + /*) A + pE (17), 

and £ is the most general solution of 

p*pg + /xV 2 ^ = ... (18). 

It can easily be verified that a solution of equation (17) is 

1 [X -f 2/x .A 

(n — - —™— I U^ — - J2j I # ... . . , , „ (It/). 

r \ p J v ; 

There will be solutions for 77, £ similar to (16), but the three solutions for £ 77, £ 

must be such that 

dP , dm , tit , 

->- 4. -1 ~f~ -^ — A . . . (20^ 

«# rty/ dz x J 

The left-hand member of (20) is, from (16), 



v ^ + ^ + % + * 





- 3 



and from (19) and (17), V 2 (f> = A. Hence (20) is satisfied if 






I -, ■— |— ■' " V' » « s « . » « , „ 1 ±j X /( 



§11. Write u for ~-^+ -™ ry + ~£ as before, and u {) for ~; : <f (J + -;' r ]o + -'^ £ . Then 
we shall verify that the solutions for w and u are 

^ zz: G£D W , 7J<q zz: 06qQ w ........ (22), 

in which a, a are functions of r, as yet unknown. 

Assuming these solutions, the value of E, calculated as explained in § 6, is 



E 



4:7TpS n f 1 



2n + 1 [r" +i 

where « a denotes (a) riS 



r' i+ 'J M+i (Kr)^r + r" r- ,,+ *J B+} («r)dr + — 1«. h • (23), 



a* 

Y 2 



164 



ME, J. H. JEAXS OX THE VIBRATIONS ANT) 



We can calculate the value of the integrals which occur in this expression, and the 
sum of the first two terms inside the curled brackets is found to be 



2n + 1 



K" 



r *J n+i (icr) 



r 



n. 



tea 



n-*~t 



% , {nay 



Hence we may write (19) in the form 



(p _ (j]xb ;/ -J- 



4:7TpS n 



where 



(2n + l)f a 



Alii 



(J£ = Or *J n+i (/cr) + DW W _ § (*a) 



« t- 9 I ^lr± ^J 



c 



X + 2/i 4?rp 
]rp 2 rfC 



9. Q 






D 



47T/} 



(271 + 1) fit 



a 



n + h 



We now have, from equation (16), 



d 



£ = T" (ffiS w ) + r 



4 



?rp 



d! 



LVJb 



2n + l)p 3 & 






+ 4 



and hence 






rf-r 



2-7 1 -f l)p" 3 ^"' 



* g- 



(25) : 



(26). 



§ 12. There are three boundary-conditions to be satisfied, expressing that the 
normal pressure and the two tangential tractions shall vanish at every point of the 
free surface. As Lamb # shows, these may be represented by three symmetrical 
equations, to be satisfied at the surface r = a, each of the type 



/£ 



\xA ~j- fJL y- (ru) + jul ( r - - — £ 



dx 



a 

\ di 



0. 



Substituting for £ and u from (25) and (26) tills becomes 



A, LS-i'/ ""i /X 



~d 

U/iAJ 



r ^ 8.) + r 'i i (ffiS.) 



\ 



dr 



_v/ \Atifj 



dx 



(ffiS„) 



+ ^ ~f>n + 'V)f dx[ar+ 



+ / x ^ (™o) + r 



# 



fen 



6^7 ' 



, 13. Write 



5 



cl 
dx 



(r"8 n ) 



r' 1 1 o), 



d 



So 






cfa 







(27). 



so that the right-hand members are solid harmonics of degrees n — 1 and — (n + 2) 
then 

* Lamb, he, ciL arrfe, p, 191. 
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r 

2n + 



^s» == o-~-j (^ — x) ; 



d 

(h. 



; f /« 8 -> = s£r {*"""" I ('**'/> - - '" i <'-/) x} 



From these identities it is clear that if the terms in (27) which do not depend on 
<f or «« are expanded in spherical harmonics, they will contain no harmonics other 
than a) and x> We therefore see that the form of £ may be assumed to-be 

i = pco + ©y • . (28), 

where |13 and © are functions of r. The value of v is 

w = (n^ - (n+ 1)©)S„ . ....... (29), 

whence 

= ^f - (n + 1) ©. (30). 



a 



§ 14. Substituting for £ in (27) and equating the coefficients of a and x, we obtain 
the two following equations which must be satisfied at r = a : — - 

x rJJ, H ,H _i_ r +]) rf / dG\ I d _ \ / d_ \ 

li In + 1 ^ 2» + 1 l ^ V ^ / V d >- I X dr ( ®V 

, 4-rrp(2n~ 2) /f\»-' 1 , ...„ d , „ , ^33 

+ -<^W W a « + 2»-n r_( } * ( r " +8 -») + '' I - » = ° • < 31 )> 

and a second equation of a similar kind, of which the first line can be obtained from 
the first line of the above by writing — (n + 1) for n, and the second line is 

1 d , , n x d<B, 

-2,Tl r "^^ < " - l)a °) + r dr -®' •■... (32). 

The expression which occurs in curled brackets in (31) can be transformed into 

2 {»'" (B+,, |(^7^ tlffi ))-^- ( " +,> |(^ +, ffi)} • • • ■ (33), 
while the corresponding expression in (32) is seen to be 



2 



o rf /'l tf 

di 



'17 v 

,(-~ : (r-^)j+(n+l)r"|(r-»ffi)} .... (34). 



From the value of <£, given by equation (24), 

;(r'^ffi) = GVr-U^^r) 4- (2w + 1) I>r"J.,^ (k«) ; 



d 

dr 



(r-"ffi)=-CKr-<" + »J B+ ,( K ,-), 
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Hence expression (33) becomes 

2 {CicVJ n _ f (kt) — nCKT~ h J n _ h (kt) + (2n + l)(n 
of which the value at r = a is 



l)Dr»-VT w _j (*<*)}, 



^ = 2^C/c 2 J^a (ica) + 2 [(2rc + 1) (n — 1) Da' 1 " 1 - n«-*Cie] J n .j (#ca). 

This is the value at r = a of the term, which occurs in curled brackets in 
equation (31). The value of the similar term in (32), namely expression (34), is seen 
to be 



Write 



2 = 2a*Cie 2 J rt+f (ica) — 2 (w + 1) a~ i C/cJ il+ . l (ica) 






» s 



. . (35). 
. . (36), 



>J = 0^ -f a>J , (/ca) 






then equations (3.1) and (32) become, at r 



- a, 



+ - v -«,, + a- ( ' 4+1) T (r' l+ -a ) + (2n + 1 r f 



w 







(38) : 



T3 + «"|, (r-^'S.) - (27i + I) (V - rf) . 

Now we have, from equation (26), 



. 



<*« = 



Write 



d(fii\ Airpn i ( \ 

dr /,._„ (In, + I.) ?r v u/ 



e 



Air pn, 

(2u + i)^ 



(2 ib + l)/> 3 — 4 777m 



then this last equation becomes 



^v j| i \ 



a a = c ( a -t- 



,//■ / 



Now, at r 



— a 



d 



a~ in+1) j;(r n+2 * ) — (n + 2)a 



'■/« 







dr 



. . . ( O «y ). 



:/ 



(t» y (7 



,. --(<*- D a 







) = ~ (u + ! ) a [} + a 



dr 



and equations (38) and (39) become 

_ , 4«irp (2>?, — 2) c / , d<8i\ , , . „* 
» + "^ ("o + ^7 J + fa + 2) * 



«^+(2^+l)fr5-») = . (40), 



clr 



di 



dr 



dr 







e « 3 « « « 



• • ( 41 ) : 



in which 7* must he put equal to a. 
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§ 15. Now £ is known to be a solution of equation (18), and hence, from equation 
(28), we may assume 



where 



W = l%r-hl„_ h (hr).. 



€ = ©<)»'-* J n+ , {hr) ; 



h? = f> l pl^ (42). 

It follows that at r = a we may write 

d<B, 



a f = m 



a ■ -- = B®, 



where 



d 



A = a y log (a KJ n _ h (ha)) . . 



da 



d 



* 9 



B = a j„ lo s ( a * J »+i M) • • 



rfo. 



» • » • • 



. . . (43), 



• ( 44 )- 



Hence, from (30), we have at r = a, 



d 



a-£ = n&$ - (n 4- 1) B® (45). 



f/? 



Lastly, we have from (28) 



d 



r/ 



£> - ^- n ~ ' .;:. (H3„) + <©r« + 2 - (r-<* +J >S„). 



^ 



r/.^ 



In order that (21) may be satisfied, we must have 



j- (fr"^) nr n ~ l - J (®r* +2 ) (rc + 1) r~ (w+2) = . 



. . (46). 



Substituting for :}$ and (5=t, we find that this is satisfied for all values of r if 



n 



p o + (n + 1) a = o 



• • 



• * • 



i.e., if we have, at r 
where 



(47), 



a. 



€-#$ ..... (48), 










I * 



« • * 



. (49). 



Equations (40) and (41) now assume the forms 






" 4xrp(2n- 2)e 

/r 



(n - (n + 1 ) 6) + (n 3 - 1 ) 4- (n 3 + 3w. -j- 2) 



+ (3n + l)A — (»+l)B0 



= . 



* • a 



. . (50) ; 



V + P [n — ra s + (rc s + 27i)0 + wA~ (3rc + 2) B0] = . . . (51), 

in which r must be put equal to a. 

The general frequency equation may be found at once by the elimination of ||, the 
values of WL and U being given by equations (36) and (37). 
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Points of Bifurcation. 

§16. The interest of the question lies in the position of the points of bifurcation ; 
to find these we must put p 2 = in the frequency equation. The reason why it was 
not possible to put jr = at an earlier stage will be understood by those who have 
read the former paper " On the Stability of a Spherical Nebula." In the present 
instance it is, perhaps, sufficient to say that putting p 3 = at an earlier stage would 
have led to an entirely misleading result. Upon putting p 2 = in equations (50) 
and (51) we find that the two brackets multiplying |3 vanish, and we therefore see 
that P must be treated as an infinite quantity of the order of 1/p 2 . 

Expanding these brackets as far as p 2 , and then putting p~ ~ 0, we find that the 
two equations become 



where 



*i - WVi = . . . . (52), x 3 + \Bp% = . . . , (5.'}), 

%r (It 

(n - 1) (2n + I) 3 , pa" f 2 CM - 1) ?>n + 1 



2/2 



2/iTrp ' /x, [ (2n + 1 ) (2-/i + 3) 2 (?t -f 1) I 

pa 3 ?i 



fi 2(2n + l)(2?i + 3) 
The equation giving points of bifurcation is, of course, 



ItJ ^ ~\ Xolj -\ ~~~ ^ • a • • • • ' • • • (04), 

The values of x } and ,r 2 are found, after some simplification, to be 
x y = - a*J n+i (*a) -f 2C j a*/cJ n _,j (k«) — ncC h K3 n ^ h (kci) — iOj^ilj:.. (a%I /l+ i (fca)) 



A* 



2 (^-1) (2 ? i+l) ^ 
+ ■-■ a ' s J w _j (^a) , . . . . . , (55), 



# 3 = - a*J n+i . (*a) + 2a-C/c 3 J, l+ , (/ca) — 2 (n + 1) a~*C/cJ n+ .». (fca) . . (56). 

1 
Now, it has already been seen that C = — — - (p. 164). If we substitute this value 

for C, write x for /ca, and simplify equations (55) and (56) as far as possible, we have 



JL i (.'V 



^ J,,„ W + *J!>±i^ j,„ 5 W _ ^=i^_tl> J-ti (x) . (57): 



~l A 4" 2^6 * , x 2 (?£ ~j-2)yyv / r \ 
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while the value of y y and y % may he written in the form 



Vi 



pa 



o r 

< 



4h s + 20?i 2 + 21» + 7_ 2 ( » - 1) (2 n + 1)» /x 1 

'" + w 2 ~ " X + 2u j ' 



(2w + I) (2» + 2) (2» + 3) 



(59), 



2/3 



/Oft" 



71 



/a 2 (2» + 1) (2« + 3) * 



(60), 



The equation giving points of bifurcation can now be found by substituting these 
values in equation (54). 

§ 17. This equation will have roots corresponding to the different integral values of 
n 9 n = 0, 1, 2 . . . ; these determine points of bifurcation such that the critical 
vibrations are of orders n = 0, 1, 2 . . . respectively. 

Of these the points of bifurcation of zero order are of no importance. The reason 
is exactly similar to that given in the case of a spherical nebula (§ 28 of the paper 
already quoted), namely, that a point of bifurcation of order n = does not indicate 
a departure from the spherical shape. We therefore will only discuss values of n 
different from zero. 

Case of /x = 0. 

§ 18. Before discussing the general form assumed by equation (54), it will be well 
to consider the simple case of /x = 0. Putting /x = 0, we obtain from equations (57) 
and (58) 



lAJ \ 






P 



Wvu + l \ X )' 



Referring to equations (52) and (53) we see that the equation giving points of 
bifurcation is 

<J//,-f* y£) ~~~ V e * • - » 4 . » ' • * (01 ). 



The lowest roots of various orders other than zero are 



n 

x 



1, 
4*49, 



2, 
5-76 






6*98, 



4, 

8*18 



9-37, &c, 



the roots continually increasing with n. Thus the first point of bifurcation is given 
by x = 4*49, and the critical vibration is of order n = 1. 



Case of /x .Different from Zero. 

§ 19. The general equation in which li is not put equal to zero is much more 
complicated than equation (61), which has just been considered. If we write \i n for 
JV + i (x)/J n ,_± (x)> it will be seen that the equation giving points of bifurcation of order 
n is of the form 



u u + l = an algebraic function of x and of (X + 2/x)//x. 
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To obtain approximate numerical solutions, my plan lias been to draw graphs of 
the functions u m and in this way obtain a graphical solution of the ecjuations for 
different values of /x. There is no difficulty in drawing graphs of the functions u fl ; 
these are trigonometrical functions, and we have 






1 



x 



cot X, 



while the successive its are connected by the relation 

2n -f 1 1 



u 



il + 1 



%l 



To save space I have suppressed all details of this somewhat tedious part of the 
work. The results for n = 1, 2, 3 are given in the following table :-— 

Lowest Values of x. 





/X=0. 


11/ "0'/v» 


/X — A. 


n=l 
n = 2 


4 • 49 
5*76 
6-98 


5*6 
6-8 


4-0 
5-4 
6-7 



For large values of n it will be found that equation (54) reduces approximately to 
x. 2 = 0, and hence that for any value of [i the lowest value of x is slightly less than 
the corresponding value in the case in which /x = (X 

jf/ie ifyVsi Pom£ q/* Bifurcation. 

§ 20. It therefore appears that the first point of bifurcation may be safely assumed 
to be of order n = 1. The value of x for which it occurs will have some value 
between 4*0 and 4*49, according to the value of /x/A. Now x = k«, and the value of 
k 2 is Airp^ftk + 2/x). Hence the first point of bifurcation is approximately given by 

4tt/5% 2 f4*00 3 = 16'00 5 when /x = X, 
H2^ [4*49 3 = 20-1(3, when /x == 0. 

In equation (3) we supposed this point of bifurcation to be given by 

y/> s a 3 Ao = <£• 

In our present analysis we have already taken y = 1 ; if we take (X + 2 m) to be 
identical with our former X , we see that the actual values of (j> are roughly 

(/> = 1'60, when /x = 0, </> = 1 8 27 ? when /x == X* # 



* It will be found that the first point of bifurcation is given, with great accuracy, by the single 
equation p% 2 /(A + lp) = 1 • 6 for all values of /x between and A, This is of interest, as showing the 
relative importance of A and /x in maintaining stability. As might be foreseen, the importance of /x 
•relatively to A increases as n increases, and for 71= oo , the factor A + f/x must be replaced by A 4- 2/x. 
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We have now found a closer approximation to the value of <j> than that which was 
given in § 3, and have obtained the additional information that instability first enters 
through a vibration of order n = 1. It must, however, be borne in mind that these 
results are only true of the special and somewhat artificial case specified in § 8. 

Comparison with the Case of a Spherical Nebula. 

§ 21. It will be seen that the general argument of § 3 will apply to the case of a 
gaseous planet or nebula if X be taken to mean the pressure in the gas. In this case, 
however, the laws of distribution of density and pressure are not independent. If 
the gas is in conductive equilibrium throughout, the planet or nebula must be 
supposed to extend to infinity, and for these conditions the criterion of stability 
was worked out in the former paper already referred to. Calling the elasticity 
of the gas k, the first point of bifurcation was found to be reached when the function 

U f —-*-— attains a certain finite value. Now U, p vanishes in comparison with p , 

the mean density, so that writing a for the radius of the nebula, and X for the mean 
pressure (X = /cp ), we have, at this first point of bifurcation 



2irp 2 a?/k 



CO . 



Comparing this with the general result obtained in § 3, we see that in this extreme 
case the value of (j> becomes infinite. This result is only of importance to the present 
investigation as showing the tendency of a concentration of density about the centre. 
It seems to show that as the density becomes more concentrated about the centre, 
the value of <[> may be expected to increase. We are therefore led to expect that in 
general <£ will have a value rather greater than that found for it upon the 'assumption 
of homogeneity of density. 

Recapitulation and Discussion of Results. 

§ 22. It will be well to recapitulate our results before attempting to draw any 
deductions from them. 

We consider a spherically symmetrical mass of solid, liquid, or gaseous matter. 
We denote the radius of this by a, the mean density by p Q , and the mean value of X 
by X , where X denotes an elastic constant or the pressure of the fluid, according as 
the matter is solid or fluid. We have seen that the stability of this dynamical system 
depends upon the value of the function yp 2 a 2 /\ , a pure number. When y = (i.e., 
when we deal with artificial matter which is totally devoid of gravitation) there can 
be no doubt that the system is stable. We have seen that a point of bifurcation 
occurs when the number y/o 3 a 2 /X has a certain value (f>. It has not been proved in 
the present paper that an exchange of stabilities accompanies this point of hi furcat ion, 

/i 4J 
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but it will be seen that, with slight alterations, the proof of the exchange of stabilities 
for the spherical nebula, which was given in § 28 of the earlier paper, can be made to 
apply to the present case. Admitting this, it appears that the spherical system 
which is at present under discussion will be stable so long as yp^(vj\ is less than <£, 
and becomes unstable so soon as yp °CL 2 /X exceeds <j>. 

§ 23. The next question is as to the exact value of <£, and as to the vibration through 
which instability enters at the point of bifurcation. To the first part of the question 
we have not been able to obtain a very definite answer. This matters the less, since 
the numerical data which would have to be used in making any applications of our 
results are not themselves very definite. On the whole, the uncertainty in the value 
of cj> is not much greater than the uncertainty in the value of the numerical data (or, 
what comes to the same thing for our present purpose, the uncertainty in our 
knowledge of the law of compressibility and distribution of density in the planets of 
our system), 

The general argument of § 8 showed that $ must, except in extreme cases, be 
comparable with unity. We then examined an artificial case : a planet in which the 
density and elasticity were constant throughout — this system being made mechanically 
possible by introducing an external field of force, of amount just sufficient to annul 
gravitation in the equilibrium configuration. For this system p was, of course, 
taken equal to p, the uniform density, and X was taken to be equal to X -f- 2/x in the 
notation of Love, or m + n in the notation of Thomson and Tait. The value of <f> 
depends, of course, on the ratio /x/X or n/m. For fi/k = we found <j> = 1 '6 ; for 
fji/X = 1 we found <j> = 1*27 ; for intermediate value of /x/X we saw that the value of <£ 
was intermediate between these two values. 

The planets to which we wish to apply our results do not possess uniform density : 
it is almost certain that in every case the mean density is much greater than the 
surface density. The general argument of S 3 shows that there is still a point of 
bifurcation corresponding to a value of </> which is comparable with unity, but affords 
no evidence as to the change which a concentration of density will effect in the value 
of (f). We therefore examined a case in which there is an infinite concentration of 
density — the case of a spherical nebula extending to infinity — -and found, that in this 
extreme case the value of <j> becomes infinite. It therefore seems probable that a 
concentration of density is attended by an increase in the value of <f>. As a working 
hypothesis we shall assume for the planets of the solar system the uniform value 
<f> = 2. It must be left to the reader to form a judgment as to the amount of error 
involved in this assumption, but it will, perhaps, be admitted that results depending 
upon it will at least be right as regards order of magnitude. It will be seen later 
that considerable variation in the value of <j6 is possible before the astronomical 
evidence which we are going to bring forward is seriously invalidated. 

§ 24. As regards the nature of the vibration through which instability of the 
spherical configuration enters, we are able to come to a more definite conclusion. In 
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each of the cases referred to in the last section this vibration is found to be one of 
order n = 1, i.e., one in which the displacement at every point is proportional to the 
first harmonic. This is the result which we should naturally expect — just as we 
expect a mass of liquid to become unstable through long surface waves sooner than 
through short ones. We shall, therefore, suppose it to be true of the planets in 
general. It is conceivable that planets could be artificially constructed for which this 
assumption would not be true, but, at present, since we have not a complete 
knowledge of the structure of the planets and are therefore compelled to make some 
assumptions, it seems as if the assumption just made is far and away the best to take 
as a working hypothesis. 



Application to the Nebula n Hypothesis. 

Theoretical Conclusions. 

§ 25. In the former paper, already referred to, the suggestion has been put forward 
that the instability of a nebula, sun or planet, which, upon the nebular hypothesis, is 
supposed ultimately to result in the ejection of a satellite, may be largely brought 
about by a gravitational tendency to instability of the kind we have been investi- 
gating. Let us, for the moment, take an extreme hypothesis, and imagine that this 
agency is the preponderating agency, and that the rotational tendency to instability 
may be disregarded in comparison. 

Upon this hypothesis let us consider the case of an approximately spherical planet 
or sun which is known to have thrown off a satellite. Before the ejection of this 
satellite commenced, the primary mass would have an approximately spherical form, 
for which p °a%/\ would be below the critical value <j>. When this critical value is 
reached, a divergence from the spherical form occurs, and a series of new processes 
begins. We are not now concerned with the details of these processes, but they 
must be supposed ultimately to result in the ejection of a satellite. It must be 
noticed that we are not supposing the primary to be devoid of rotation — for this 
would be inconsistent with the ejection of a satellite — but are supposing the rotation 
to be so small that the rotational tendency to instability is small in comparison with 
the gravitational. 

If we suppose one or more satellites to have been ejected, and the primary to have 
regained an approximately spherical form, the new value of p 2 a 3 /\ must be less 
than cfj. Now every satellite of which we have any knowledge, in our own system 
at any rate, is small in comparison with the primary. A legitimate inference seems 
to be that the ejection of a satellite is only a small part of the life-history of the 
primary. We shall not, however, need to make any assumption so definite as this, 
but shall suppose only that the values of p 0P a, X for the primary after ejection are 
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nearly equal (i.e., equal, except for a multiplying factor of, say, 1^ or 1^) to the 
values of the same quantities before the process of ejection commenced. 

All this may be summed up, with sufficient accuracy for the present investigation, 
in the statement that when the ejection of a satellite is taking place the value of 
/> 3 a 2 /X must be nearly equal to <£. 



Co?iclu$ions Tested by the Solar System, 

§ 26, To a certain extent the solar system supplies material for testing this 
conclusion. Let p , a denote the present mean density and radius of any member of 
the system which is known to have thrown off a satellite, and let X be the mean 
value of X, whether this denotes a pressure or an elastic constant, then we may 

write 

/) 2 a 2 /X = 0<j>. 

At the time at which the satellite was ejected the same equation ought to hold 
with 6 nearly equal to unity. Unfortunately we have no knowledge as to the 
changes which will have taken place in p , a, X since the ejection of the last satellite. 
We shall, for a rough approximation, neglect these changes altogether, take 0=1, 
<f) = 2, and examine to what extent the equation 

holds for the solar system as it now stands. 

It will be remembered that there are three sources of error in this equation : — - 

(i.) We are neglecting the effect of rotation in bringing about the ejection of a 
satellite. 

(ii.) We are neglecting the changes which have taken place since the ejection of 
the last satellite. 

(iii.) We are using an arbitrarily chosen value for <j>, and applying this to every 
planet, while we know, from the difference in the physical constitutions of the planets, 
that the value of </> must be different for each. 

These three sources of error would each be serious if we were attempting to get 
accurate results, but as our calculations are necessarily only of the roughest kind, we 
may be content to neglect them. 

§ 27. In the accompanying table the masses and radii of the sun, and of those 
planets which possess satellites, are given in the first two columns ; the value of 
mass/ (radius) 2 is given in the third. The units are so chosen that the earth is 
measured by unity in each of these columns. Venus is included for the sake of 
comparison, although the existence of a satellite is extremely doubtful, 
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Observed. 




Calculated upon the hypothesis 

of the present paper 

(0 = 2). 


(i) 

Mass. 


(2) 
Eadius. 


(3) 

Mass 
(Radius) 2 


(4) 

Coefficient A . 

Unit = 10 11 absolute 

= 10 s grammes weight per 

sq. centim. 


Venus .... 
Earth .... 
Mars . . . 
Jupiter .... 
Saturn .... 
Uranus .... 
Neptune . . . 


315,000 

0*8 
1-0 

o-i 

300-0 
90-0 
14-0 
16-0 


109 
1-0 
1-0 
0-5 
11-0 
9-0 
4-0 
4-4 


26 
0*9 
1-0 
0-4 
2-5 
1-1 
0-9 
0-8 


2700 
3-2 
4-0 
•6 
25-0 
5-0 
3-2 
2-6 



If our hypotheses give a fair account of the facts the numbers in this third column 

will be proportional to \/\^>. Assuming for <£ the uniform value <j) = 2, Ave can 
calculate the actual values of X , and these are given in the fourth column. 

§ 28. Knowing nothing about the variation in X , we shall be content as a 
preliminary hypothesis to suppose it to have the same value for each, planet. 
Combining this with the hypotheses already formulated, we notice that \Z\tff> ought to 
have the same value for each planet, as therefore ought also the function 
mass/ (radius) 2 , which is tabulated in column (3). 

It will be seen at once that there is a certain amount of uniformity about the 
numbers in this column, but it requires some consideration to determine how much 
significance is to be attached to this uniformity. 

Now, apart from any hypothesis as to how the solar system originated or reached 
its present configuration — i.e., regarding the solar system as a fortuitous collection of 
bodies of varying sizes— Ave should expect the mean density to be greatest in the 
greatest planets. We should, therefore, expect the quantity (mass) /(radius) 3 to be 
more variable than the radius. In other words, we should, a priori, expect less 
uniformity in the third column than in the second. Judged by this criterion, the 
uniformity of the numbers in the third column would be very significant. Further, 
the variation in these numbers is of the kind we should expect. For instance, it is 
known that the density of Jupiter is very much greater near the centre than near the 
surface ; we should accordingly expect a large value of <£, and therefore a large entry 
in the third column. The same argument would apply to the Sun, but the physical 
constitution of the Bun is probably so different from that of the planets that there 
could be no surprise at the Sun figuring as an, exceptional case. Another exception 
is that of Mars. Part of the discrepancy might, perhaps, be attributed to the 
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smallness of the planet, but the figure in the fourth column would seem to suggest 
that rotational instability must have played a large part in the creation of the 
Martian satellites. 

If, on the other hand, we begin by regarding the planets not as a fortuitous 
collection of bodies, but as a series of satellites all ejected from the same primary, the 
case is different. For here we should expect the smaller planets to have cooled more 
than the heavier ones, and therefore to be at a lower temperature. Against this 
must be set the fact that the heavier planets will probably have the greatest 
concentration of density about the centre, and the greatest mean pressure. The first 
consideration tends to increase the value which we should expect for the mean 
density of the smaller planets as compared with that of the greater ones; the second 
consideration tends in the opposite direction. We can hardly profess to estimate the 
relative weights of these two considerations with any approach to accuracy ; perhaps 
it is best to revert to the argument given in the last paragraph, while bearing in 
mind that the approximate equality of our numbers- may become considerably less 
significant as soon as the question of relative temperature is taken into account. , 

§29. We now consider the evidence afforded by the absolute value of our figures, 
After allowing for the exceptional cases, it appears that the value of X for the earth 
and for most of the planets is about 4 X 10 11 . In other words, if we suppose these 
planets suddenly to resume the molten state, while retaining their present mass and 
radius, the spherical form will be stable or unstable according as the mean value of 
X + 2/x is greater or less than 4 X 10 u . In the molten state we may take jx = 0, 
and the value X = 4 X 10 11 corresponds to a value equal to about half of that of 
steel, for which X = 8 '3 X 10 11 . If, however, we attempt to trace the history of a 
planet backward in time, we cannot suppose the mass and radius kept constant : the 
mass may be constant, but the radius will increase. Under these conditions we find 
that the critical value of X will be inversely proportional, to the fourth power of the 
radius, and will, therefore, be somewhat less than the value )i=4X 10 u . It would 
be extremely difficult to form a reliable estimate of what this corrected value for X 
ought to be, and equally difficult to estimate what would be the actual value of X for 
molten material similar to that of which our planets must have been composed when 
in the molten state. Our argument is that the two values of X are at least of the 
same order of magnitude, and probably equal, except for inaccuracies in our calcu- 
lations. 



Comparison of the Rotational and GravitaMonal Hypotheses, 

§ 30. We may conclude this part of our work by comparing two extreme 
hypotheses : the first referring the phenomena of planetary evolution solely to 
rotational, and the second solely to gravitational instability. 

Given the approximate values of X and p for a planet, and the fact that this 
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planet has thrown off a satellite, the former hypothesis leads to a certain inference as 
to the angular momentum of the system ; the latter to an inference as to the radius 
of the primary. The former hypothesis leads to no inference at all as to the size of 
planets which are to be expected — they are as likely to be of the size of billiard balls 
as of the size of the planets of our system— while the latter leads to no inference as 
to the angular momentum of the system, but presupposes it to be small The 
contention of the present paper is that the inferences which are drawn from the 
former hypothesis are not borne out by observations on the planets of our system, 
while those which are drawn from the latter are borne out as closely as could be 
expected. The true hypothesis must of necessity lie somewhere between the two 
extremes which w r e are comparing, and our evidence seems to show that it is much 
nearer to the latter (gravitational) than to the former (rotational). # 



Stresses and Vibrations in the Earth. 

§ 31. It has already been seen that in dealing with a gravitating sphere of the size 
of the earth it is necessary to take into account terms which are omitted by Lord 
Kelvin and others — the terms which introduce into our equations the gravitational 
tendency to instability. 

It is of some importance to know whether the existing solution for the vibrations 
and displacements of the earth would be altered to an appreciable extent by the 
inclusion of these terms. The general frequency-equation which is given on p. 167 is 
too complicated for manipulation, and is, moreover, open to the objection that it does 
not represent the facts of the case ; for, inside the earth, the strains caused by 
permanent gravitation cannot legitimately be treated as small, f 

§ 32. Considerations of a general nature will, however, give us some insight into 
the question. In an imaginary earth, in which X, jjl are infinitely great, the 
gravitational terms will be of no importance in comparison with those representing 
the elastic stresses. The true solution will, therefore, become identical with the 
classical solution in which the gravitational terms are neglected. For smaller values 
of A, jjl the error will become appreciable, and if X, /x continue to decrease this error 
will become infinite as soon as the first point of bifurcation is reached ; for at a point 
of bifurcation the application of an infinitesimally small external force will produce a 
finite displacement in the solid. For intermediate values of X, jx the error will be 
small if X, //, are great compared with the critical values of X, /x at the point of 
bifurcation, and great if X, /x are near to these critical values. 



* 



In addition to the inference as to the size of the planets, the hypothesis of gravitational instability 
leads to a further inference as to the distances of the fixed stars. This has been discussed in my former 
paper, " On the Stability of a Spherical Nebula " (§ 48). and here also the results seem to agree with 
observation as closely as could reasonably be expected. 

t Chree, ' Camb. Phil. Soc. Proc.,' vol. 14, or Love, ' Elasticity,' vol. 1, p. 220. 
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§ 33. The- most reliable evidence as to the actual values of X, fi is to be obtained 
from the phenomena of earthquake propagation. # From the "time curves" given in 
the British Association Report presented at the 1902 meeting, there seems to be 
little doubt that the so-called "large-waves" are propagated merely through a thin 
crust on the earth's surface, while the "preliminary tremor" is propagated in a 
sensibly straight line through the earth itself. The average velocity of propagation 
is found to be about 9 8 7 kiloms. per second, and this is independent of the length of 
the path. The inference is that (X ~f- 2/x)/p is nearly constant throughout the earth's 
interior, and that its value is about (9*7 X 10 5 ) 2 ' or 9*4 X 10 1L . If we suppose the 
mean value of p to be 5\5, this gives for the mean value of X + 2//,, 

\ + 2jjl = 51 7 X 10 11 . 

Now, the critical mean value of A.+.2/X, which corresponds to the first point of 
bifurcation has already been seen to be about 4X 10 u . It would, therefore, appear 
that the error introduced in the classical solution for the displacements and stresses is 
appreciable, although not great — it is probably comparable with the error to which 
attention lias already been attracted by Oh'KEE. t 

FlGUKE OF THE EaIITH. 

Theoretical Conclusions. 

§ 34. From the evidence of the last section it will be seen that there is an over- 
whelming probability that the values of the elastic constants of the earth are such 
that a state of spherical symmetry would be one of stable equilibrium. 

Whether or not the earth is at present in a state of spherical symmetry is a 
different question ; various indications and, in particular, the inequality in the 
distribution of land between the two hemispheres of the globe suggest that it is 
not so. 

Now. even if the material of the earth is at the present moment of sufficient 
strength to maintain a spherical configuration in spite of the gravitational tendency 
to instability, it does not seem, probable that it has always been so. Looking back- 
wards in time we must come to a stage in which the material of the earth was plastic, 
and, further back still, fluid. At this time the value of X would be much smaller 
than its present value, and, as already pointed out in § 29, would probably be about 
equal to the critical value for the planet at that period of its existence. There would, 
therefore, seem to be a sufficient reason for considering the possibility that the earth, 
at the moment at which consolidation set in, was not in a state of spherical symmetry. 
Let us examine some of the consequences of this conjecture, 

* Professor Milne has kindly assisted me in this question, 
t Log, eit, ante. 
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It is easy to see that enormous stresses would be set up in the interior of the earth 
after consolidation. An equilibrium configuration depends in general upon the 
compressibility of the material, and a configuration which was one of equilibrium for 
the compressibility which obtained at the moment of solidification would not remain 
so after the incompressibility and rigidity of the material had increased by cooling. 
If we suppose the earth to cool in an unsymmetrical configuration the stresses set up 
will soon become very great. In fact, Professor Darwin has shown that the stresses 
which would be produced by the weights of our continents in an earth initially 
homogeneous (i.e., by an. irregularity of less than a thousandth part of the radius) 
would be so great that the material would be near the breaking point. # 

We must therefore suppose that as the earth cools and the elastic constants change 
there will be a series of ruptures resulting from the stresses set up in the interior. 
The configuration. will become approximately spherical (spheroidal if rotation is taken 
into account) as soon, as the point of bifurcation is passed. 

The fact that the ultimate configuration is readied only as the result of a long 
succession, of ruptures puts the whole question outside the range of exact mathe- 
matical treatment. We can, however, see that the filial configuration (disregarding 
rotation) will probably not be quite spherical, but will retain traces of the initial 
unsy mm etrical con figu rati on . 

§ 35. Before we can attempt to decide whether or not the earth shows traces of a 
process such as that just described, it will be necessary to form some idea of the 
unsymmetrical configuration with which the process must have begun. We cannot 
accurately calculate the " linear series" of unsymmetrical configurations except in the 
immediate neighbourhood of the point of bifurcation. Near to this point the 
configuration is spherical except for terms proportional to the first harmonic. The 
free surface will, therefore, be strictly spherical, and it will, of course, be an equi- 
potential. but its centre will not coincide with, the centres of other surfaces of equal 
potential. If we suppose a fluid mass of this kind to solidify, and then to shrink by 
cooling, the shrinking being accompanied by a series of ruptures of the kind already 
explained, we can easily imagine that the free surface would retain an approximately 
spherical form, but that when the final state is reached this surface would not be 
quite an equipotential, and the centre of gravity would not quite coincide with the 
centre of figure. If water is placed on the surface of a planet of this kind, it will 
form a circular sea, of which the centre will be on the axis of harmonics, while the 
dry land will form a spherical cap. 

Evidence from the Distribution of Seas and Land. 

§ 36. Now this is not observed on the earth, and it could not be expected, since we 
have ignored all the agencies which have contributed to the figure of the earth, 

* 'Phil Trans./ vol. 173, 1882, p. 187. 
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except the one with which this paper is specially concerned. The question is not 
whether we observe the state just described, but whether we can detect any approach 
to this state, and this, I believe, can be done. Professor Darwin writes* : — 

"It is well known that the earth maybe divided into two hemispheres, one of 
which consists almost entirely of land and the other of sea. If the south of England 
be taken as the pole of a hemisphere, it will be found that almost the whole of the 
land, excepting Australia, lies in that hemisphere, whilst the antipodal hemisphere 
consists almost entirely of sea. This proves that the centre of gravity of the earth's 
mass is more remote from England than the centre of figure of the solid globe. A 
deformation of this kind is expressed by a surface harmonic of the first order." 

§ 37. We can carry our calculations a step further. The divergence from the 
initial configuration is only represented by a first harmonic so long as squares of this 
divergence may be neglected. If these squares are taken into account, we must 
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include a term proportional to the second harmonic as well as that proportional to 
the first. This process of successive approximation might be continued to any extent, 
so that a complete series of unsymmetrical configurations might be calculated in the 
manner explained in my former paper, f We may, however, be content to stop at 
the second harmonic . The free surface will now be of the form 

* Q. H. Darwin, ' Phil. Trans./ vol. 173, 1882, p. 230. 
f 'Phil Trans.,' A, vol, 199, p. 41. 
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r = a + fti-Pi + %?2 * (63), 

and we therefore examine whether any traces of the second harmonic term can be 
found in the earth's surface. Now, if we take a 2 positive in equation (63), the 
equation is that of the pear-shaped curve which was found on p. 46 of this earlier 
paper. This differs from the spherical shape mainly in possessing a protuberance— 
the stalk end of the pear — of which the centre is on the axis of harmonics. Traces 
of this protuberance may, I think, perhaps be found in the Australian continent, the 
arrangement being that shown in fig. 1. It is true that the centre of Australia does 
not coincide with the antipodes of England, but the discrepancy becomes less when 
we take into account the enormous region of ocean shallows which lies to the east of 
Australia. 

# The discrepancy can be further reduced by taking the rotation of the earth into 
account. When the rotation of the earth was greater than at present the ellipticity 
of the earth's surface would be greater, and the transition from this to the present 
ellipticity would take place through a series of ruptures similar to those already 
described. The rotation (assumed small) of the pear can be allowed for by adding a 
term — /3P ; 2 to the right-hand side of equation (63), this representing a second 
harmonic deformation having the axis of rotation for axis of harmonics. 

The present rotation of the earth can similarly be represented by a term — /3T' 2 , 
where y S / < j3. The equation to the present surface of the sea may accordingly be 
taken to be 

r = «' ~ PP' S , 

and hence the height above the present sea-level of the surface of the primaeval 
rotating pear, if restored, would be 

K - a' ) + a v Pi + « 2 P 2 ~ 08 - P) P' 2 . 

It will be found that the effect of the rotational term (/J — ft) P ; 3 is to move the 
theoretically predicted Australia nearer to the equator of the earth, and to change its 
shape from a spherical cap to a sphero-conic] 

Again, we should expect the highest land to be on the axis of harmonics, and, 
therefore, in or near England. Here, again, the agreement of facts with theory 
might be closer if we could suppose the continent, which geology shows to have 
existed at one time in mid- Atlantic, to be restored to its former position. But the 
agreement of facts with theory can only be expected to be of the roughest kind, and 
we must always bear in mind that our theory does not lead us to expect that the 
present figure of the earth will be pear-shaped, but only that it will resemble a pear 
disfigured by a long series of ruptures. 

* Added January. 3, 1908, I am indebted to the referee for suggesting this addition, 
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Evidence from the Distribution of Earthquake Centres, 

§ 38. It can be seen that the earthquake regions of the world have a reference, as 
regards their distribution on the earth's surface, to this pear-shaped figure, and this, 
again, must be considered as evidence. 

Let us first examine the facts. Milne divides the earthquake-areas of the globe 
into twelve distinct regions, and a map of these is given in the c British Association 
Report' for 1902.^ These regions are given in the following table. The first figure 
denotes the number of large earthquakes which have occurred in these regions in the 
three years 1899-1901. The earthquakes from the three regions printed in italics 
were small in comparison with the others. In the last column is given the 
approximate latitude of the centre of each region, referred to Greenwich as pole (the 
latitude of Greenwich being taken to be + 90°). 

Table of Earthquake Regions. 
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Now, it will be at once noticed that for most of these regions the latitude is small. 
If we weight the regions according to the corresponding number of earthquakes, 
giving half- weight to the small earthquakes in regions H, 1, J, we find as the mean 
of the numerical values of the latitude about 20° ? whereas if the regions were 
distributed at random we should expect the mean latitude to be (\tt — 1) radians, or 
about 33°. We therefore see that the earthquake regions tend to lie near the 
equator of our pear. The evidence can be put in a more striking way as follows : — 
Exactly half of the surface of the globe is of a latitude less than 30°. The half for 
which the latitude is less than 30°, measured from Greenwich as pole, was responsible 
for 156 earthquakes; the remaining half was responsible only for 42, of which 28 
were the small earthquakes from regions H, I, J. There is, therefore, no doubt that 
the principal earthquakes tend to emanate from points near to the equator of the 
supposed pear. 

Now, if we look back to fig. 1, we see that this is equivalent to saying that earth- 
quakes occur where the il slope ?; in the figure of the earth is steepest. This conclusion 
is the same as that to which the British Association Committee were led from a 

* * Brit. Assoc, 72nd Report/ Belfast, 1902, " Seismological Investigations," p. 4. 
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consideration of the actual figure of the earth, and it is that which might naturally 
be expected. The theory put forward in this paper may, perhaps, suggest a reason 
why these regions should lie approximately on a great circle of the earth, and why 
this great circle should approximately divide the earth into two hemispheres of sea 
and land. 



Summary and Conclusion. 

§39. In conclusion it may be well to summarise those parts of the paper which 
refer to the figure of the earth. 

We saw that at the moment of solidification the earth might be either spherical 
(except in so far as it was deformed by its rotation) or pear-shaped. Our theoretical 
calculations and our knowledge of the constants of the earth at the time of solidi- 
fication were not sufficiently accurate to enable us to decide which of the two 
alternatives is the more probable. The shape of the earth, whether spherical or pear- 
shaped, could not be maintained long against the enormous strains which would be 
set up in the earth as the process of cooling proceeded, and this shape would gradually 
give place to an approximately spherical shape, the change in shape being produced 
by a long succession of ruptures. The suggestion of this paper is that the earth, in 
spite of this series of ruptures, still shows traces of a pear-shaped configuration. 
Such a configuration would possess a single axis of symmetry, and this, it is suggested, 
is an axis which meets the earth's surface somewhere in the neighbourhood of 
England (or, possibly, some hundreds of miles to the S.W. of England). Starting 
from England we have in the first place a hemisphere which is practically all land ; 
this would be the blunt end of our pear. Bounding this hemisphere we have a oreat 
circle of which England is the pole, and it is over this circle that earthquakes and 
volcanoes are of most frequent occurrence. If we suppose our pear contracting to a 
spherical shape we notice that it would probably be in the neighbourhood of its 
equator that the change in curvature and the relative displacements would be 
greatest, and hence we should expect to find earthquakes and volcanoes in Greatest 
numbers near to this circle. Passing still further from England we come to a oreat 
region of deep seas— the Pacific Ocean, the South Atlantic Ocean, and the Indian 
Ocean : these may mark the place where the "waist" of the pear occurred. Lastly 
we come, almost at the antipodes of England, to the Australian continent and the 
shallow seas which lie to the east of it ; these may be the remains of the stalk-end of 
the pear. 

§ 40. It may, I am afraid, be thought that the hypotheses upon which the paper is 
based are too speculative and the results, consequently, too uncertain. In defence it 
may be said that the object of the paper is not so much to establish new doctrines as 
to point out possibilities, and that these possibilities seem to be of such a kind that it 
may be useful to keep them in mind in discussing questions connected with the figure 
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and structure of the earth, as well as the more general questions of planetary 
evolution. 

In conclusion I have to express mv indebtedness to Professor G. H. Darwin and 
Professor A, E H. Love for advice and assistance which 1 have received from them. 

[Note,— Added February 20th, 1908. While the above paper was in the press, 
Professor W. J. Sollas read a paper before the Geological Society in which the 
Figure of the Earth was discussed from a geological standpoint. Professor Sollas 
had arrived, from an. examination of the geological features of the earth, at a 
conclusion very similar to that to which X had been led from theoretical 
considerations : he had detected an axis of symmetry, other than the axis of rotation, 
in the earth's figure, and expressed the opinion that "the pear-shaped form, now that 
it was pointed out, became obvious to mere inspection i it was a geographical fact, 
and not a speculation." 

The axis of Professor Sollas' pear does not, however, coincide with that which I 
tentatively put forward in the above paper, and the object of this note is to accept 
the alteration suggested by his paper. The conclusion reached in his paper is that 
the axis of symmetry of the pear-shaped figure passes through a point of latitude 
and longitude about 6° N. by 30° E. Thus Africa— the continent whose mean height 
above sea-level is greatest— must be taken to be the centre of the " Land 
Hemisphere " in fig. 1 of my paper, while the protuberance which formed the stalk 
of the pear is submerged in the Pacific Ocean, which, now forms the " Water 
Hemisphere/' Almost the only remaining evidence of the existence of this 
protuberance is the fact that the axis of the pear coincides with the earth's greatest 
diameter. The great circle of earthquake-centres suggested in § 38 of my paper is 
to be replaced by the line of Pacific folding ; this approximately forms a small circle 
(of radius about 80°) which almost coincides with the proposed great-circle in the 
northern hemisphere. Further details of Professor Sollas' view will be found in his 
paper ("The Figure of the Earth," 'Quart. Journ. Geol. Soc,,' vol. lix., Part 2). 

The fact that Africa is surrounded by a belt of seas, and tin's again by a belt of 
land before the Pacific is reached, points perhaps to a bodily subsidence of the blunt 
end of the pear, the circle of fracture having possibly been the line of Pacific folding. 
Such a fracture would, of course, displace the centre of gravity of the pear, and 
probably this would account not only for the feature just mentioned, but also for the 
non-appearance of the protuberance. It will be noticed that the smallness of the 
latitude of the extremities of the axis (6°) agrees well with the theory of planetary 
evolution put forward in §§ 25-30 of the present paper. J 
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